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AN  EFFICIENT  COMPUTATIONAL  PROCEDURE  FOR  SOLVING  THE 
MULTI-FACIUTY  RECTIUNEAR  FACILITIES  LOCATION  PROBLEM^ 

Henrik  Juel  and  Robert  Love 

INTRODUCTION 

A problem  that  occurs  frequently  In  the  area  of  facilities  location 
concerns  the  optimal  location  of  several  new  facilities  In  relation  to  a set 
of  existing  facilities.  The  problem  that  is  addressed  here  Is  the  one  in 
vhich  distances  can  be  considered  to  be  rectangular.  A typical  situation 
occurs  when  one  or  more  new  machines  must  be  located  on  a plant  floor  on 
which  there  are  existing  machines  and  storage  areas.  The  floor  is  arranged 
in  a rectangular  grid  such  that  material  movement  must  take  place  parallel  to 
two  orthogonal  axes.  The  objective  is  to  minimize  the  total  material  flow 
costs  between  the  new  machines  and  the  existing  machines  and  storage  areas. 
Other  applications  include  the  location  of  service  facilities  Inside  office 
buildings,  the  location  of  factories,  warehouses  and  postal  stations  in  urban 
areas  and  the  design  of  piping  and  wiring  circuits.  Further  discussics  of 
appllc.itlons  are  found  in  reference  9. 

A practical  application  involving  the  location  of  two  new  facilities  in 
relation  to  three  existing  manufacturing  departments  is  given  by  Love  and 
Terex. ^ For  this  relatively  nmall  problem  it  was  possible  to  use  linear  pro- 
gramming as  the  solution  technique  using  a formulation  developed  by  Wesolowsky 
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and  Love  and  Cabot  et.  al. 

^Sponsored  in  part  by  the  United  States  Army  under  Contract  No.  DAAG29-7  5- 
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Center  under  these  grants  was  made  possible  through  support,  in  part,  from  the 
National  Science  Foundation,  other  U.S.  government  agencies  and  the 
Wisconsin  Alumni  Research  Foundation. 

*Visltlng  Professor  in  the  Department  of  Management  Sciences,  University 
' of  Waterloo,  Waterloo,  Ontario  in  1976-77. 


For  larger  problems,  however,  the  linear  progrannlng  formulation  of  the 

I 

location  problem  may  give  rise  to  extremely  large  numbers  of  variables  and 
equations.  As  an  example,  consider  a problem  with  m-350  existing  facilities 
and  n«70  new  facilities.  The  linear  programming  primal  has  2mn  + n^  “ 53,900 
variables  and  mn  + n(n-l)/2  - 26,915  constraints.  Even  the  more  efficient 
linear  programming  dual  has  26,915  bounded  variables  and  70  constraints.  The 
direct  search  approach  of  Pritsker  and  Ghare^  is  more  appealing  from  a compu- 
tational point  of  view,  since  much  larger  problems  can  be  solved  than  with  the 
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linear  programming  formulation.  However,  Rao  and  Wesolowsky  have  provided 
counter-examples  proving  that  the  Pritsker  and  Ghare  procedure  nay  not  reach 
optimality  in  all  cases.  The  present  work  may  be  viewed  as  an  optimizing 
implementation  of  their  approach. 

In  this  paper  we  state  and  prove  necessary  and  sufficient  conditions  for 

optimality  in  all  cases.  (These  conditions  were  indicated  by  Rao  using  the 

theory  of  linear  programming.)  A solution  algorithm  is  formulated  incorporating 

the  conditions,  and  computation  results  are  given.  These  indicate  that  the 

algorithm  is  very  rapid  when  large  numbers  of  new  facilities  are  not  located 

at  a single  point.  We  also  give  a new  lower  bound  on  the  optimal  solution  to 

■ultl-faclllty  problems.  This  bound  is  generalized  to  account  for  t distances 
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which  are  discussed  by  Love,  and  Love  and  Morris. 

PROBLEM  AND  NOTATION 

The  multi-facility  rectilinear  facilities  location  problem  In  one  dimension 
can  be  formulated  as: 

minimize  f(x)  - [ I - a | + j £ ''211 '*1  " *1  • * 

J-1  1-1  ^ J ^ 1-1  J-i+1  ^ J 


-I- 


where  x ~ (x,,  x.,  x },  d ia  the  number  of  existing  facilities,  and  n Is 

i « n 

the  number  of  new  facilities,  “nj  non-negative  weight  between  existing 

facility  1 and  new  facility  J,  for  l£l£">>  ''21j  non-negative 

weight  between  new  facilities  1 end  J,  for  1 f.  1 < J £ n.  a^^  is  the  (one- 

dioensional)  location  of  existing  facility  1,  for  1 £ 1 £ m.  Xj  Is  the  (one- 

dlmenslonal)  location  of  new  facility  j,  to  be  determined,  for  1 £ J £ n.  It 

fulfices  to  consider  the  problem  In  one  dimension,  since  an  unconstrained  prob- 

q 

lem  In  more  dimensions  can  be  decomposed  into  one-dimensional  problems 

It  Is  convenient  to  assume  a,  < a,  < ...  < a and  to  define  w.,,  = w,, . for 

12  m 2J1  21J 

1 £ 1 < J £ n.  Let  us  denote  by  a current  solution  a situation  In  which  each 
new  facility  coincides  with  some  existing  facility.  It  suffices  to  consider 


current  solutions,  since  an  optimal  solution  exists  which  is  a current  solution. 
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For  a current  solution  x = (x, , x-,  ...,  x ),  let  s,  denote  the  index-set 

1 z n K 

of  new  facilities  currently  coinciding  with  existing  facility  k,  that  is: 

Sj^  - for  1 £ k £ m.  Finally,  let  f (x;y)  denote  the  direc- 

tional derivative  of  the  convex  function  f at  the  point  x = (xj^,  x^,  ....  x^) 

, ^ , I ^ ^ M r\,  ^ -,±  f(x+hv)  - t (x) 

In  the  direction  y “ (y,,  y^,  ....  y ),  that  Is,  f (x;y)  * lim  'T . 

^ ^ " hHH  ^ 

For  the  objective  function  of  the  multi-facility  rectilinear  facilities 

2 

location  problem,  at  a current  solution  x,  we  obtain; 


f(x;y)  - I 
k-1 


I b y + I I + I “ikjlyjl 

jcs  J J jcs.  its  jts 


k k 
l<j 


where  b^  denotes  the  force  exerted  on  new  facility  J by  the  facilities  not 
coinciding  with  It,  that  Is, 


bj  " bj(*)  " 1 Wjjj»fgo(*j~aj)  + I w^jj  sign(xj-x^)  for  jcsj^,  l£k£m. 
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since  the  objective  function  la  piece-wise  linear,  the  directional  derivative 
is  Just  the  increase  In  the  value  of  the  objective  function  caused  by  a suf- 
ficiently small  perturbation  of  the  current  solution.  That  Is, 
f'(x;y)  “ f(x  + y)  - f(x),  provided  that  the  components  of  y are  small  enough 
In  Mgnltude  that  the  restriction  of  f to  the  line  segment  connecting  x and 
X + y Is  linear. 

NECESSARY  AND  SUFFICIENT  CONDITIONS  FOR  OPTIMUM 
Since  the  objective  function  f Is  convex,  a necessary  and  sufficient  con- 
dition for  a current  solution  x to  be  an  optimal  solution  Is  that  the  direc- 
tional derivative  f’(x;y)  be  non-negative  for  all  directions  y.  This  condition 
Is  not  particularly  operational,  and  to  transform  It  into  a useful  result,  we 
need  the  following: 


Pact: 


for  all  y-(yj^,y2 y^^) 


‘ ;"k  ’ ' i"k  ^ ' 

A^  of  a : I I b I < [ I w + J w 
JtA^  U\ 


if  and  only  if,  for  all  subsets 
means  IcSj^-A^). 


Proof: 

The  only-lf-part  follows  by  setting  Pj  " * Pj  “ 0 for 

For  the  If-part  we  need  2 Intermediate  results.  First  suppose  t ■ 

■1®  ^ 1*1  c be  any  real  number  satisfying 

t + C £ 0 and  * + c £ y^  for  1*1  ll>*  vector  y'  be  given  by  y^  >■  y^  + c 

for  JtAj^  and  yj  “ for  J^Aj^.  Then  by  direct  computation: 
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Next  suppose  i - max  {y^cjcsj^)  >0.  - ^{JcSj^ryj  - r).  Let  c be  any  real 

number  satisfying  x-e  ^ 0 and  x-c  ^ y^  for  j^A^-  Let  the  vector  y'  be  given  by 
7j  " Xj  “ c for  JcAj^  and  yj  » y^  for  Then  by  direct  computation 

Now,  for  any  given  y we  can  apply  these  results  repeatedly  with  c > 0,  together 
with  the  assumption  of  the  If-part,  obtaining  gj^(y)  “ gj^^(O)  + non-negative  terms  ^ 0. 
This  fact,  in  conjunction  with  the  previously  given  expression  for  the  directional 
derivative,  proves  the  following. 


Result: 

A current  solution  x is  an  optimal  solution  if  and  only  if  for  all  subsets 

,,  for  1 < k < m. 

1 ■ ' ' — . zn  XK.  r 

JcAj^  J JcAj^  leSj^( 


A.  of  8.  (x):  I j;  b (x)|  1 I I ''211'^  ^ 

' * JcA,,  J JcA^  les^(x)-A^  J j£. 


'Ikj' 


For  use  in  the  algorithm,  2 special  cases  of  the  result  are  given. 
1.  If  the  subset  “ (j),  then  the  condition  is: 


l»>z(x)I  1 I 


+ W, 


lcSj^(x) 

i+J 

If  this  condition  is  not  satisfied,  an  improved  current  solution  is  obtained 
by  moving  the  single  new  facility  J from  a^^  to  “^-sign  (bj  (x) ) ' 

2.  If  the  subset  Aj^  “ Sj^(x),  then  the  condition  is: 

I I ^(x)l  1 I ''ikj- 

JeS|^(x)  ^ JcSj^(x)  ^ ^ 

If  this  condition  is  not  satisfied,  an  improved  current  solution  is  obtained 
by  moving  the  whole  cluster  of  new  facilities  in  Sj^(x)  from  a^  to 


•b-signf  I bx(x)V 

IJcs^(x)  ^ J 
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ALGORITHM 


On  the  basis  of  the  preceding  section,  the  following  algorithm  is  pro- 
posed: 

1.  Find  a current  solution. 

2.  Improve  the  current  solution  by  moving  new  facilities  one  at  a time, 
until  no  further  Improvements  are  possible  (using  special  case  1 of  Result). 

3.  Improve  the  current  solution  by  moving  whole  clusters  of  coinciding  new 
^*cllltles,  one  whole  cluster  at  a time,  until  no  further  Improvements  are  pos- 
sible (using  special  case  2 of  Result). 

A.  Check  all  subsets  of  each  cluster  of  coinciding  new  facilities  (using 
Result).  If  moving  some  subset  Improves  the  current  solution,  make  that  move 
and  restart  step  A.  Otherwise,  an  optimal  solution  has  been  found. 

This  algorithm  has  been  Implemented  much  like  the  one  by  Pritsker  and  Chare.  Thus, 
the  initial  current  solution  is  found  by  solving  n single-facility  problems,  ig- 
noring the  weights  between  pairs  of  new  facilities.  Any  move  consists  of  moving 
one  or  more  new  facilities  from  the  current  location  to  the  location  of  a neigh- 
boring existing  facility.  After  each  move  the  sre  updated,  rather  than 

recomputed  from  the  basic  data. 

COilPUTATlONAL  RESULTS 

The  algorithm  was  programmed  in  Fortran,  and  some  test  problems  were  run  on 
a Unlvac  1110  computer.  The  computation  times  reported  here  represent  CPU  time 
exclusive  of  compilation  and  reading  or  generation  of  input  data.  The  2 test  prob- 
lems given  by  Pritsker  and  Chare  took  .05  seconds  each.  Eight  test  problems  with  m>50 
and  n”10  were  generated,  the  weights  being  random  numbers  between  0 and  99. 

These  problems  took  an  average  of  .9  seconds  each.  Most  of  the  time  was  spent 
In  step  A,  since  the  typical  optimal  solution  involved  most  or  all  of  the  10 


new  facilities  forming  a cluster  on  or  near  >25*  ^ problems,  In  6 cases 

the  optimal  solution  was  found  at  step  2,  In  1 case  at  step  3 and  in  1 case  at 
step  For  the  typical  problem,  step  4 just  performed  the  optimality  test  as 
intended. 

Test  problems  with  ni“350  and  n“70  were  generated  with  biased  weights  to 
avoid  the  excessive  and  probably  atypical  clustering  of  new  facilities.  The 
weights  were  random  numbers  between  0 and  9,  except  that  for  1 j £ 70 

and  5(j-l)  + 1 £ 1 Jl  5j  were  random  numbers  between  0 and  W-1.  Fourteen  prob- 
lems with  W=300  took  on  the  average  1.2  seconds  each.  In  one  case  the  optimal 

solution  was  found  at  step  2,  In  4 cases  at  step  3,  and  In  9 cases  at  step  4. 
Fourteen  problems  with  W-400  took  on  the  average  1.1  seconds  each.  In  no  cases 
was  the  optimal  solution  found  at  step  2,  In  9 cases  at  step  3 and  in  5 cases 
at  step  4.  In  these  28  problems,  most  of  the  time  was  spent  In  step  2.  In 

some  problems  with  a smaller  value  of  W,  an  excessive  amount  of  time  was  spent 

in  step  4,  and  the  program  was  terminated  before  completion.  For  example,  in 
a problem  with  W*200,  the  solution  after  step  3 involved  a cluster  of  21  new 
facilities.  Probably  more  than  20  minutes  would  have  been  spent  just  to  check 
all  the  subsets  of  this  cluster.  (The  computer  can  check  roughly  1000  subsets 
per  second.) 
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APPENDIX  - PROPERTIES  OF  THE  I FUNCTION 

P 


Let  F (x)  denote  the  objective  function  of  a multi-facility  location  problem 
P 


with  1 distances  In  d dimensions: 
P 

n ■ 


n-1  n 


%(*)  - I I -in  + I I “21J  Vi  ~ 

^ J-1  1-1  P J ^ 1-1  j-i+l  P ^ J 


where  x—  (x^*  *2*  •••*  *n^  * ^J  " ^^J 1 * ^J  2 * * * * * ^J  d^  fur 

•1  ■ ^*11*  “12 *ld^  for  1 < 1 < m and  t^Cz)  - Zj)  ■ 


I 

k-1 


1/p 


Let  hj^  denote  the  minimal  value  of  the  corresponding  kth  one-dimensional 
subproblem: 


n m n-1  n 


Then  a lower  bound  on  the  optimal  value  of  the  d-dimensional  problem  is  given 
by  the  t^-norm  of  the  vector  h »■  (hj^,  hji  -..t  : min{Fp(x))  ^ Tl.is 

generalizes  a result  by  Prltsker  and  Chare  (6),  who  consider  the  case  p=2,  d=2. 
The  result  follows  from  the  triangle  Inequality  for  t^-norms. 


Fd  f n m 

III  I 

lk-lU-1  1-1 


1/p 


1 I I' 

J-1  1-1 


llj 


^''lljl*jk'“lkl 

n-1  n 

+ Z I ''21l'*ik‘ 

1-1  J-1+1  ^ 

•*JkI 

P 

1/p 

' d 

n-1  n 

[ d 

1 

^ I’‘jk’“lkl‘* 
lk-1  ■* 

+ I 1 «2iJ 

1-1  J-1+1  ^ 

JJ^k-^Jkl 

1/p 


- Fp(x)  . 
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